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Abstract. We investigate the evolution and the statistical properties of the Ly a absorbers of the intergalactic medium (IGM) in 
the largely unexplored redshift range z = 0.5 - 1 .9. We use high-resolution (R > 30 000) UV (STIS) and optical ( VLT/UVES and 
Keck/HIRES) spectra of nine bright quasars with z cm < 1.94. The Ly a lines detected in the lines of sight (LOS) towards these 
quasars are evaluated with a software package which determines simultaneously the quasar continuum and the line profiles. 
The main results for the combined Ly a line sample are summarized as follows: 

1. The evolution of the number density of the absorbers can be described by the power law $j oc (1 + z) y ■ The number density of the low 
column density lines (A'hi = (10 1290 - 10 1400 ) cirT 2 ) decreases with decreasing z with y = 0.74 ± 0.31 in the interval z = 0.7 - 1.9. 
A comparison with results at higher redshifts shows that it is decelerated in the explored redshift range and turns into a flat evolution for 
z — > 0. The stronger absorbers (A'hi > 10 1364 cirT 2 ) thin out faster (y = 1.50 ± 0.45). The break in their evolution predicted for z = 1.5 - 1.7 
cannot be seen down to z = 0.7. On the other hand, a comparison with values from the literature for the local number density gives a hint 
that this break occurs at lower redshift. 

2. The distribution of the column densities of the absorbers is complete down to A'hi = 10 12 90 cirT 2 . It can be approximated by a single power 
law with the exponent /? = 1.60 ± 0.03 over almost three orders of magnitude. /? is redshift independent. 

3. The Ly a lines with lower column densities as well as the higher column density lines show marginal clustering with a 2cr significance over 
short distances (Av < 200 kms~' and Av < 100 kms~', respectively). We do not see any difference in the clustering with either column 
density or redshift. 



The distribution of the Doppler parameters has a mean value of b ■■ 
not change significantly with z. 



(34 ± 22) km s . This value is typical for the analyzed region. It does 
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1. Introduction 

It is traditional to decompose spectra of the Ly a forest into in- 
dividual absorbers that can be fit with Voigt profiles, each with 
a redshift z, a column density of neutral hydrogen Nh i and a 
Doppler parameter b. Examining a sufficiently large number 
of LOS to background quasi-stellar objects (QSOs), one can 
derive the distributions of these parameters as well as depen- 
dences between them and in this way deduce characteristics of 
the Ly a forest. Observers have often focused on the evolution 
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of the number density of Ly a lines, usually approximated by 
the power law 
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(Sarg ent et all fl980). At high redshift (z > 1.5), the evolu- 
tion of the strong lines (log Nhi > 14.00) is steep (y = 
2 - 3), transiting into a fiat evolution (y = 0.1 - 0.3) 
at lower redshift. This abr upt slow-down in the evolution , 
first detected by obser vers ( Bahc all et al. L fl99ll llmpev et all 



ll996tlWevmann et a 111993) and later also seen in simulations 
jTheuns etalH l998: iDave et allll999l) . is suspected to happen 



atz= 1.5-1.7. 

However, since the apparent break in the evolution rate 
seems to occur just at the transition from high-resolution opti- 
cal spectra to low-resolution spectra taken with the HST, the lo- 
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QSO 




m B [mag] 


PG 1634+706 


1.34 


14.9 


PKS 0232-04 


1.44 


16.6 


PG 1630+377 


1.48 


16.5 


PG 0117+213 


1.50 


16.1 


HE 0515-4414 


1.73 


15.0 


HE 0141-3932 


1.80 


16.2 


HE 2225-2258 


1.89 


16.3 


HS 0747+4259 


1.90 


15.8 


HE 0429-4901 


1.94 


16.2 



Table 1. Parameters of the quasars 



cation of this break needs to be confirmed using high-resolution 
UV spectra. In a previous study, using UV observations of the 
extremely bright z em = 1.73 QSO HE 0515-4414 with a res- 
olution of 30000, it appeared that the slow- down in the evo- 
lution rate did not occur earlier than z ~ 1 ( I Janknecht et all 
120021) . However, at low redshift there are few absorbers, and 
more LOS are necessary. 

In this paper we use UV and optical high-resolution data 
of nine bright quasars to explore the Ly a forest in the red- 
shift range z = 0.5 - 1.9. We believe that this is the first study 
of the lower column density absorbers in the redshift region 
z < 1 .5 using spectra with R > 30 000. The paper is organized 
as follows: in Section 2, we explain our selection criteria for 
the quasars, describe the observations and define the redshift 
ranges that we investigate. Section 3 illustrates how the Ly a 
lines in the spectra were detected and modelled. After mention- 
ing some interesting individual features in the individual LOS 
in Section 4, we discuss our results in Section 5. The main re- 
sults are summarized in Section 6. 

2. Observations 

2.1 . Selection of the quasars 

Table^gives an overview of the quasars whose spectra we se- 
lected to analyze the Ly a forest. We list their emission red- 
shifts z em and their apparent B magnitudes mj. We selected 
QSOs with emission redshifts suitable to show the evolution 
of the Ly a forest at 1 < z < 2. In Table|5]we give the spectral 
resolution and the approximate S /N of the spectral regions that 
we examined. The resolution of the spectra of all nine quasars 
is high enough to fully resolve all Ly a lines. The S/N varies 
from one spectrum to the next, and within each spectrum, so 
the minimum observable column density depends strongly on 
the spectrum and the wavelength. 

2.2. Spectra and data reduction 

Because the Ly a line is shifted into the optical for z > 1.5, a 
complete analysis of the Ly a forest at median and high red- 
shifts requires both optical and UV spectra. For two of the nine 
quasars, HE 0515-4414 and HS 0747+4259 (discovery of both 
objects published in Reimers et al., 1998; for the Ly a forest 
in the LOS to HE 0515-4414 alone see Janknecht et al., 2002), 
data in both spectral regions have been obtained: in the UV 



with HST/STIS with the E230M echelle mode (R ~ 30000), 
and in the optical with UVES (R ~ 50000) for HE 0515-4414 
and with Keck/HIRES (R ~ 38 000) for HS 0747+4259. 

The sample was extended with the UVES spectra of three 
bright quasars from the Hamburg ESO Survey taken in the ser- 
vice mode wit h the VLT. These thr ee were HE 0141-3932, 
HE 2225-2258 JWisotzki et all l2000l) and HE 0429-4901 (not 
yet published). The last was taken with a slit width of d — 0.8", 
the other two with d - 1.0". The Ly a lines detected in these 
spectra strengthen the statistics for 1.5 < z < 2.0. 

Finally, the QSO sample was supplemented in the UV 
(0.5 < z < 1.5) with data from the STIS archive. For 
this purpose the spect ra of four quas ars — PKS 0232-04 
dShimmins et all Il966l). PG 1630+377 dNoguchi et all E98O) 
and PG 01 17+213 Tschmidt & Greerl Il983l) from an observ- 
ing program of B. Jannuzi (HST propos al 8673/cycle 9) and 
PG 1634+706 dSchmidt & GreenL Il983l) from two programs 
of S. Buries (7292/7) and B. Jannuzi (8312/8) — were down- 
loaded from the STIS data base. These data were also taken 
with the E230M echelle mode resulting in the same resolution 
applicable to the other STIS spectra. 

Table|2]lists the observing dates, exposure times and quality 
informations for all the QSO spectra. 

The data reduction was performed using the pipelines of 
UVES and STIS on their respective data. The Keck data of 
HS 0747+4259 was reduced by J. O'Meara using an internal 
pipeline and the data reduction software package of T. Barlow. 
For the HE 0515-4414 STIS spectra another rework had to be 
made: the radiative flux was in places smaller than zero because 
an inaccurate correction of the noise background of the CCD 
chips was used in the pipeline. The problem was solved by 
measuring the signal needed to make the minimum flux in satu- 
rated lines zero. A slightly reduced background value (varying 
from order to order, and extrapolated to orders that lacked sat- 
urated lines) was then subtracted from the raw data flux, result- 
ing in positive flux values for the entire spectrum, apart from 
random photon noise. 

The vacuum and the baryocentric corrections to the wave- 
length scales were performed with the data analysis package 
MIDAS. Multiple exposures of each QSO were added, weight- 
ing each flux value F^j by the inverse variance (noise) crj 2 . 

2.3. Definition of redshift ranges 

The sensitivity variations with wavelength of the individual de- 
tectors slightly reduce the redshift regions in which the Ly a 
forest can be examined. For HS 0747+4259, we defined a min- 
imum signal to noise per pixel, averaged over each spectral 
order, for inclusion in the analysis. This is because the sen- 
sitivity of the HIRES chip decreases strongly with decreasing 
wavelength. Choosing (S /N) > 2, we ignored the first three or- 
ders for HS 0747+4259. Analogously, the first two orders for 
HE 0429-4901 were left out. The spectrum of PG 1634+706 
contains a Lyman limit system, making any detection of Ly a 
lines below A < 1865 A impossible. We will discuss below 
how we separate Ly a from other Lyman lines. 
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QSO 


Telescope / spectrograph 


Observing date 


Exposure time [s] 


R 


S IN per pixel 


PG 1634+706 


HST/STIS 


May 5, June 26, 1999 


26400-29 000 


30000 


5-50 


PKS 0232-04 


HST/STIS 


Feb. 6, 8, 2001, Jan. 19, 2002 


41900 


30000 


4-15 


PG 1630+377 


HST/STIS 


Feb. 27, Oct. 8, 12, 2001 


34 100 


30000 


5-11 


PG 0117+213 


HST/STIS 


Dec. 31, 2000-Jan. 12,2001 


42 000 


30000 


4-12 


HE 0515-4414 


HST/STIS 


Jan. 31-Feb. 2, 2000 


31500 


30000 


~ 10 




VLT/UVES 


Oct. 7, 2000-Jan. 3, 2001 


31500 


50000 


10-50 


HE 0141-3932 


VLT/UVES 


July 19, Aug. 14-24, 2001 


39 600 


40000 


~ 25 


HE 2225-2258 


VLT/UVES 


June 17-28, July 14, 15, 2001 


41 800 


40000 


~ 25 


HS 0747+4259 


HST/STIS 


Sep. 6, 12-16, 2001 


54200 


30000 


3-10 




Keck I/HIRES 


Feb. 28, Mar 1, 2001 


5 400 


38 000 


3-12 


HE 0429-4901 


VLT/UVES 


Feb. l,Mar 18, 19, 2001 


10 800 


50000 


~ 8 



Table 2. Log of the observations used in this study 



QSO 


AA [A] 


z range 


Az 


AX 


n 


PG 1634+706 


1865 -2790 


0.534- 1.295 


0.761 


1.665 


195 


PKS 0232-04 


2280 - 2941 


0.876- 1.419 


0.544 


1.309 


128 


PG 1630+377 


2279 - 2980 


0.875 - 1.451 


0.577 


1.396 


118 


PG 01 17+213 


2279 - 3009 


0.875 - 1.475 


0.600 


1.460 


160 


HE 0515-4414 


2278 - 3260 


0.874- 1.682 


0.808 


2.034 


220 


HE 0141-3932 


3061 - 3384 


1.518 - 1.784 


0.266 


0.745 


97 


HE 2225-2258 


3057 - 3478 


1.515 - 1.861 


0.346 


0.979 


130 


HS 0747+4259 


2140 - 2970 


0.760- 1.443 


0.683 


1.615 


110 




3115 -3484 


1.562- 1.866 


0.304 


0.864 


79 


HE 0429-4901 


3188 - 3538 


1.622- 1.910 


0.288 


0.830 


88 


I 






5.176 


12.897 


1325 



Table 3. Evaluated redshift ranges and absorption distances of the quasar LOS and numbers of Ly a lines n detected in the 
indicated spectral regions 



We also excluded redshift regions within Av ~ 2 500 - 
5 000 kms -1 of the emission redshifts where the ionization of 
hydrogen is enhanced by the proximity effect. We calculated 
the maximum redshifts, given in Table [3] from the absolute 
magnitudes of the QSOs. 

For HS 0747+4259, there is a gap of approximately 145 A 
between the UV and the optical spectra, while for HE 0515- 
4414 there is an overlap of roughly 50 A. In this overlap, we 
used only lines which can be identified in both spectra (nine 
altogether) and we chose the line parameters from UVES be- 
cause they are more accurate. Six lines are only found in one 
spectrum (three in the STIS and three in the UVES spectrum, 
respectively), all with column densities < 10 13 - 20 crrT 2 . We 
suppose that these are predominantly artefacts. 

Absorbers with constant comoving number density and 
constant proper sizes will have a changing density per unit 
redshift. We define the absorption distance AX so th at non- 
evolv ing absorbers have a constant density per unit X (Tvtler, 
Em. Assuming the ACDM model, 

1 + z 

AX = Az (2) 



jMisawaet all [2002). where £2m and Qa are the density pa- 
rameters for matter and for dark energy, respectively. 

Given a set of QSO spectra, in general there are multiple 
LOS which sample a specific redshift range. We take this into 



account by converting the z into X values and summing the AX 
values to obtain the total absorption distance. 

Table [3] shows the wavelength and redshift regions 
(columns two and three) that we investigated for all quasars 
after removal of the proximity zones and low S /N orders. The 
fourth and fifth columns give the computed redshift path Az 
and the absorption distance AX. The total absorption distance 
for all nine examined quasar LOS is AX = 12.897. In the sixth 
column we give the number of detected Ly a absorption lines, 
1325 altogether. 

In Fig.^ we show the redshift intervals of the LOS, while 
in Fig. |3 we show the distribution of all 1325 Ly a lines as 
a function of redshift and column density, giving qualitative 
information on the contents of the sample. It can be clearly seen 
that the redshift interval analysed in this paper is well sampled. 



3. Voigt profile fits 

The analysis of absorption lines detected in quasar spectra is 
usually performed in three steps. First the lines have to be iden- 
tified as certain atomic transitions that have been redshifted to 
the observed wavelengths. Then they are modelled. Finally they 
are accepted or rejected for the line sample. 

For the line identification we used the empirical fact that 
the Hi Doppler parameter b is rarely smaller than 15 kms 1 
and never < 10 kms -1 , while for metal lines almost always 
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PG 1634+706 
PKS 0232-04 
PG 1630+377 
PG 01 17+213 
HE 0515-4414 
HE 0141-3932 
HE 2225-2258 
HS 0747+4259 
HE 0429-4901 



0.5 1.0 1.5 2.0 



Using the redshifts of the strongest Hi systems, we could 
search for metal lines. We found many additional metal lines at 
the redshifts of the common Civ 1548/1550 doublets. 

The absorption line list derived in this way still contains 
a few unidentified lines, noise dips interpreted as absorption 
lines by mistake, unidentified blends, and wrong identifica- 
tions. However, at these low redshifts, and especially in the 
high-resolution UVES spectra with their high S /N we guess 
that such errors are less than 5% of the total list. Besides, 
the unidentified or wrongly identified lines are typically those 
with very low column densities and/or low Doppler parameters, 
which are either rejected by applying the significance level (see 
below) or do not belong in our Ly a line sample. On the other 
hand, the blends are a problem that is easily underestimated. 

We fit the Ly a lines with two different programs: the lines 
found in the HST/STIS and VLT/UVES spectra of HE 0515- 
Fig. 2. 4414 were fit with th e MID AS software package FITLYMAN 

Distribution of all 1325 Ly a lines as a function of redshift and col- dFontana & Ballesterl 1 19951) . which needs a normalized spec- 
umn density trum as input. In order to obtain this, the effective background 

continuum of the quasar was defined with MIDAS by specify- 
ing points in line-free regions and fitting these by a polynomial. 
b < 10 kms -1 . In most cases, it can a priori be discriminated An improvement to FITLYMAN is the program CANDALF 
between Hi and metal lines by inspecting the line width. developed by R. Baade. In CANDALF, the continuum is deter- 

The identification strategy was as follows. After rejecting mined simultaneously with the line fitting procedure. With this 
the interstellar lines, the strongest Ly a absorbers (i.e., the lines parallel approach, the continuum can also be defined reliably 
with the highest column densities) were picked out. Having in spectral regions where it is hidden by a high line density, 
identified a Ly a line at the wavelength A Lya with a Hi column CANDALF was used for the analysis of all quasar spectra ex- 
density log Nm ^ 14.00 (the exact lower bound depending on cept for HE 0515-4414 whose spectra were already fit before 
the S IN), the associated Ly j5 can in principle be found at the the development of CANDALF. 

same redshift. For higher column density systems, further tran- Extensive investigations of complex absorption line ensem- 
sitions of the Lyman series (with monotonically decreasing line bles, several with partly blended atomic transitions and lines of 
strength) can be detected. multiplets fit simultaneously, were carried out with both pro- 



Fig.l. 

Visualization of the analyzed redshift ranges of our quasars 
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grams to check their consistency. The comparison of the ionic 
parameters determined with FITLYMAN and with CANDALF 
showed that they were consistent within 1 cr in almost all cases. 
CANDALF converges more easily compared with FITLYMAN 
in regions with high line densities, finds the global fit minimum 
with larger probability, and its fit errors are a little bit smaller; 
however the choice of the fitting program does not have any 
significant influence on the results of the analysis of the whole 
Ly a sample (Section^. 

FITLYMAN and CANDALF both assume Voigt profiles 
convolved with the instrumental profile. The programs adjust 
three independent parameters per line to minimize the x 1 '- 
z, Nhi and b which comprises the thermal and the turbulent 
broadening of the lines. The general fit strategy, especially for 
blends, was to start with a single line and to add further compo- 
nents as long as the^ 2 decreased significantly. Different atomic 
transitions of the same ion were fit with identical values for 
Nm, b and z. This is particularly important for the saturated 
Ly a lines (log Nm £ 13.50 - 14.00) lying in the flat part of 
the curve of growth where log Nni is very inaccurate. In addi- 
tion the number of subcomponents of an absorber is harder to 
determine at higher column densities. However, the fragmen- 
tation of an absorber into several weak absorbers can be seen 
in Ly /3 and other Lyman series lines because these lines are 
usually not saturated and give more accurate parameters. 

Since the lower column density lines can hardly be distin- 
guished from the photon noise, it is usual practice to define the 
significance level 



cr w 



(3) 



where W is the observed equivalent width of the line and crw 
is the lcr error of W ( Youn g et all Il979i) . Usually the ap- 
proximation <Tw ~ with (S IN) being the average sig- 
nal to noise in the region of the consi dered absorption line is 
used JTvtler et all Il987t ICaulel Il989h . This describes the er- 
ror caused by the determination of the co ntinuu m. However it 
neglects the fit error (Semb ach & Sava ge. 1992). Here we per- 
formed an extended estimation of SL that includes the fit error: 



A 



cr w 



\R(S/N) 



+ crL f . 

Fit 



(4) 



Because of this extended error estimation, we chose a com- 
paratively low selection threshold SL > 1. With this selection 
criterium about 5% of all presumed Ly a lines were rejected 
from the sample. However, we did not reject Ly a lines that 
are at redshifts where metal lines can be seen. 

The full list of all 1325 detected and modelled Ly a lines 
with their fit parameters and errors are available in electronic 
form at the Centre de Donnees astronomiques de Strasbourg 
(CDS) data base via anonymous ftp to cdsweb.u-strasbg.fr. 

4. Individual LOS 

Before we analyze the Ly a line sample from all nine quasar 
LOS in Section|5J we briefly mention some interesting charac- 
teristics in the individual LOS. 



4.1. HE 0515-4414 

In the LOS to the quasar HE 05 15-4414 a damped Ly a (DLA) 
system can be seen at z ~ 1.15. Because of this absorber 
and the brightness of the quasar (mb = 15.0 mag), it was 
already a topic of several studies with different goals by the 
Hamburg , group jde la Varga et al,[ |2QQ0[ iReimers e^ all 20.0J 
20031 IJanknecht et a ll l2002t lOuast et all l2002i l2004 12006 



Levshakov et al. , 2003 



)03) 



4.2. HS 0747+4259 

The metal line sp ectrum of HS 0747+4 259 has been studied in 
a separate paper jReimers et all E006) which concentrated on 
the Ovi absorbing systems in this LOS. 



4.3. HE 0141-3932 

HE0141-3932is an unusual quasar in several respects: its spec- 
trum contains — atypical for quasars — an extremely weak 
Ly a emission line. In addition the emission lines of differ- 
ent ions give rather different redshifts, depending on the de- 
gree of ionization (e.g., Tytler & Fan, 1992). Furthermore, the 
quasar has several associated absorption systems presumably 
built by gas ejected from the QSO. The chemical composition 
of this gas is rather atypical. These spe cial characteristics of 
HE 01 41-3932 are discussed i n detail in lReimers et all d2005l) 
and in Levshak ov et all 12005). 



5. Analysis and Discussion 

5.1. Column density distribution 

The differential distribution function of the hydrogen column 
densities f(Nni) is usually defined as the number n of Ly a 
absorption lines per column density interval A^Vh i and per ab- 
sorption distance AX (AX from (|2j; Tytler, 1987). Since our 
data contains more than one LOS, we need the sum of the indi- 
vidual AXj over all LOS: 



f(Nm) 



AN m E/ =1 AXi 



(5) 



In Fig. the distribution function for the lines of all nine 
LOS is plotted against the column density. Note that the er- 
ror bars on log / increase with increasing column density be- 
cause of the decreasing number of lines. Fig.[3]shows that the 
data points at the highest column densities have a low statisti- 
cal significance. The fits plotted in this diagram were derived 
by weighting the data points with 1 / cr 2 og ,. The absorption lines 
were binned in intervals AlogA^Hi = 0.1. As expected, the lin- 
earity of log / breaks off at low column densities. This presum- 
ably reflects a selection effect: the weak lines are hidden by the 
the noise of the spectra. 

In general, the column density distribution function can be 
well approximated by the power law 



f(N m )=A-N } 



Hi • 



(6) 
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log Mii 


Alog/V H , 


log A 


P 


12.90-15.70 


0.1 


9.4 ±0.3 


1.60 ±0.03 


12.90-17.20 


0.1 


9.3 ±0.3 


1.59 ±0.02 


13.00-15.50 


0.5 


10.0 ± 0.4 


1.64 ±0.03 


13.00-16.50 


0.5 


9.9 ±0.3 


1.63 ±0.02 



do:o points □ 
best fit for 13. CO < log N < 16.50 



Table 4. Fit parameters of the column density distribution for 
different boundary conditions. Given are (from left to right): 
the log A^h i range of the data points considered in the fit; the 
interval width; the fit parameters log A and p with their lcr 
errors. 



id mffijjj 

l 4M S 



data points □ 

best fit 'or 12.90 < log N < 15.70 

best fit 'or 12.90 < log N < 17.20 



OH 



Fig. 3. Distribution of the Ly a column densities in intervals 
A log Ajhi = 0.1, statistical error of log / and fits to the dis- 
tribution for different fit regions. For the calculation of the fit 
parameters, the data points were weighted with their error. 



where f3 is the (negative) slope and log A is the intercept of the 
logarithmic representation of f(Nni)- 

We set a completeness limit log A,hi = 12.90 for the ana- 
lyzed line ensemble, while we chose log Nhi = 15.70 for the 
upper limit (see Fig. [3}. Using these boundaries and fitting 
the model (j6jl to the column density distribution, we derived 
log A = 9.4 ± 0.3 and p = 1.60 ± 0.03. This result is consis- 
tent with other studies of the column density distribution at 
lower log Ajhi at comparable z: Dobr zvcki et al . (2002) found 
B < 1. 6 -1.7 performing a curve of growth analysis. iHu et alJ 
ill 995b determine d y8 = 1.46 for log N m = 12.30 - 14.50, and 
iKim et alJ ( 120011) derived an exponent p — 1 .70 - 1 .74 with 
varying column density regions. 

Occasionally it has been claimed that at least two 
power laws are necessary to accurately describe the 

column density distribution o ver the observed ran ge of 

column density ( Cars well et al., 1987; Giallongo et al., 1993; 
Meiksin & Madau. Tl993t IPetitiean et all Il993t IPenton et all 
2004). f we extend our fit region to higher column densities 
(log Nm = 12.90 - 17.20), the distribution does not flatten 
(J3 — 1.59 + 0.02). However, our sample is too small at large 
column densities to see minor changes in this Nm region. 

Table 0] summarizes the results for the fit parameters log A 
and p for different column density regions and interval widths. 



Fig. 4. Distribution of the Ly a column densities in intervals 
AlogAjHi = 0.5, statistical error of log / and fit to the distri- 
bution. The data points were weighted with the error for the 
fit. 



Of course the choice of the size of the intervals A log A,h i 
in which the lines are collected is somewhat arbitrary. 
Fig. |4] is based on a larger step size AlogAjHi = 0.5. For 
log Ajhi = 13.00 - 16.50, the distribution can be well approxi- 
mated by p — 1.63 + 0.02 (again weighted with of og y)- This is 
within lcr consistent with the slope for the interval width 0. 1 in 
the comparable range log Nm = 12.90 - 17.20 (only two lines 
lie in the interval log Nm = 16.50 - 17.20), which makes the 
point that the choice of the interval width does not influence 
the result. 

Frequently a dependence of the index B on z has 
been postulated by ob servers (Kimetal., 199 71 12001M2 002: 
Dave & Trinn bOOlt iHeap et all 120021 iMisawai l2002h. 
Dave & Trip J ( b00 ll) found p = 2.04 + 0.23 and IHeap et al " 
(|2002|) derived = 2.02 + 0.21 for z ~ each, while lKim etal 
COO ill calculated p = 1.72 ± 0.16 for (z ) = 1.6 and p = 
1 .38+0.08 for (z) = 2 . 1 .heifer et al J J2002ll got B = 1.41 +0.05 
for (z) ~ 2.3 and lHu et alJ J 19951) and lKim et all dl997l) found 
P — 1 .46 and p — 1 .4, respectively, both for (z) ~ 3 (all val- 



ues calculated for the low column density region log Nhi < 14 
where the number of lin es make a reliable judgement possible). 
Simulations of Theuns et al. ( 1998J), e.g., suggested a clear re- 
lationship P(z), too. These results are contradicted by the com- 
paratively low gradients of Penton et al. (2000, 2004) for the 
local universe (fi = 1.72 + 0.06 and p = 1.65 + 0.07, respec- 
tively). However, the difference might be caused by the analy- 
sis method since Penton et al. (2000, 2004) assumed a constant 
Doppler parameter for the line fits. 

In order to search for a possible dependence /3(z), we ex- 
amined the column density distribution function separately for 
three z intervals of the same width (0.5-1.0, 1.0-1.5, 1.5-2.0). 
For this purpose, the LOS to the quasars had to be divided ac- 
cording to the chosen z intervals, and the absorption distance 
AX had to be computed for the individual intervals following 
(|2j- We selected A log Ajhi =0.1 for the interval width. The dis- 
tributions for all three redshift regions and the fits to them are 
presented in Fig. |5] In each case the fits take into account the 
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Fig. 5. Distributions of the Ly a column densities in intervals 
AlogAm = 0.1, statistical errors of log / and fits to the dis- 
tributions for different redshift intervals. Only the data points 
within the range log Am = 12.90 - 15.70 are considered for 
the fits; they are weighted with the error of log /. 



points in the interval log Am = 12.90 - 15.70. Obviously, the 
fit parameters given in Table [5] are independent of the redshift 
range. The increase of /3 from z = 1.5-2.0 (J3 = 1.55 + 0.04) 
to z = 1.0 — 1.5 (J3 — 1.58 + 0.04) represents at best a marginal 
trend which is not statistically significant. 



z 


logA'Hi 


log A 


P 


0.5 - 1.0 
1.0-1.5 
1.5 - 2.0 


12.90-15.70 
12.90-15.70 
12.90-15.70 


9.1 ±0.7 
9.1 ±0.5 
8.7 ±0.6 


1.57 ±0.05 

1.58 ±0.04 
1.55 ±0.04 



Table 5. Fit parameters of the column density distribution for 
different redshift ranges. Indicated are (from left to right): the z 
interval; the log Ah i range of the data points considered in the 
fit; the fit parameters log A and /? with their lcr errors. 



z 


n 


b [kms -1 ] 


6 median [kmS -1 ] 


£>o- [kms '] 


0.5 - 2.0 


1325 


34 ±22 


28 


22.7 


0.5 - 1.0 


270 


35 ±23 


29 


18.9 


1.0- 1.5 


595 


33 ±22 


28 


23.5 


1.5-2.0 


460 


34 ±22 


28 


23.5 



Table 6. Parameters of different b distributions. We give the 
number of Ly a lines detected, the average value, the median 
of b, and on the right a parameter that represents the mode. 



5.2. Distribution of Doppler parameters 

To investigate the distribution of the Doppler parameters, we 
counted the number of lines per Doppler parameter interval in 
interval widths Ab = 5 kms -1 . Fig. |6] shows how b is dis- 
tributed across the complete line sample (without the 24 broad- 
est lines with b > 100 kms -1 ). The distribution has a typical 
form: nearly Gaussian with a maximum at b — (20-25) km s _1 , 
an average value b = (34 + 22) kms -1 lying above the maxi- 
mum, and a long tail to higher Doppler widths. A large frac- 
tion of the lines with b > 100 kms"" 1 might represent unre- 
solved blends of several components, whereas the steep de- 
cline at small Doppler parameters is caused by the S L limit 
(lines with small b values have small equivalent widths and are 
more likely to be missed) and by our Ly a selection threshold 
b > 10 kms -1 . 

The average and median of b fit the result s from other stud- 
ies in comparable redshift regions very well: Kim_gt^lJ (12001) 
found ^median = 28 kms -1 for(z) = 1.6, while Kim et al. (2002) 
derived b = 32.6 kms -1 for (z) = 2.2. 

Frequently a dependence b(z) is claimed in the sense of a 
Doppler widt h increasing with decreasing redshift. In the high 
redshift range Lu et al. A 19961) found b = (23 + 8) kms 1 for 



(z) = 3.7. iKim et alJ Jl997) derived b median = 26 kms 1 for 
(z) = 3.35 and medi an = 30 kms -1 for (z) = 2.31, respectively, 
and lHuet al.ldl995h found b = (28 + 10) kms -1 for <z> = 2.9. 
The average value of this study for (z) = 1.31 appears to sup- 
port this t rend. Going to z = 0, the evolution of b possibly 
continues: Pent on et alJ J2OOO) derived b = (38 ± 16) kms -1 
for the local universe. 

Therefore we investigated whether an evolution of b with z 
can be determined within the z range examined here. In Fig.0 
the b distribution is shown for three different z ranges of the 
same width, while Table|6]gives parameters of the distributions. 
The highest average value b - (35 + 23) km s -1 for the interval 
z - 0.5 - 1.0 lies only marginally above those in the higher 
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Fig. 6. Distribution of the Ly a Doppler parameters in inter- 
vals Ab = 5 kms -1 . 

redshift regions. We see no evolution in the b mean and median 
over redshifts 0.5 to 2. In general, the questi on of whether the 
Doppler width changes with z is not resolved: Ki m et alJ (2002) 
also did not see any evol ution of b in the inter val (z) = 3.3 — * 
2.1 in their data, nor did lDave & Tripdd200lh (b = 25 kms 1 
for (z) - 0.2), whose results do not fit the evolution pattern 
outlined above. 

In addition to the mean and median, we also list a parame- 
ter that represents the mode, or most common b value, because 
this tends to be less sensitive to the S /N and the details of line 
and continuum fitti ng. We list values for t he paramete r brr = 
1 .0 574 ■ b^v fromlHui & Rutledgg] dl999l).lTvtler et alJ (p004 ) 
and LTena et all d2Q05t> found that the iHui & Rutledge ( 1999) 
fitting formula with a single parameter b a = (23.6±1.5)kms 
gives excellent fits to the b distribution at 1.5 < z < 2.4. Here 
we obtain essentially the same value for our two higher redshift 
bins, 1.0 - 1.5 and 1.5 — 2.0, and a smaller value at lower red- 
shifts 0.5 - 1.0. The b value distribution is significantly wider 
than the fitting function for the two lower redshift bins, and 
hence these two b, T values are not well determined. Both these 
two lower redshift bins show many more lines with b < 15 
kms' 1 than we see a t z > 1.5. We know, for example from 
lKirkmanetaII(l2005l) . that as the S /N of a spectrum drops, as 
it does here dramatically at low z, we will measure an excess of 
both low and high b values. Hence, all the changes that we see 
may arise from the changing S IN. We do not claim to detect 
any change in the intrinsic b distribution with z. 

5.3. Clustering of Ly a absorbers 

The clustering properties of the baryonic matter can be investi- 
gated by computing the correlation between the Ly a absorbers 
along the LOS. The degree of clustering of the Hi absorbers 
can be expressed by the two-point velocity correlation function 



f(Av) 



«obs (Av) 
n sim (Av) 



1 



(7) 



(Sarg ent et ail fl980) where « D b s , the number of observed Ly a 
line pairs in a given velocity separation bin Av, is compared 
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Fig. 7. Distribution of the Ly a Doppler parameters in inter- 
vals Ab — 5 kms -1 , varying redshift ranges. 



with the number of expected line pairs n s i m derived in the same 
velocity difference bin in a spectrum with randomly placed 
lines. For two absorbers at the redshifts z\ and zi, the veloc- 
ity interval at mean redshift in the rest frame is given by 



Av 



c(Z2 - Zl) 
1 T 2 



(8) 



We determined n s ; m from Monte Carlo simulations by dis- 
tributing our line sample — in accordance with the derived 
number density evolution (see Subsection I5.4i — randomly 
over the spectrum, counting the line pairs for various velocity 
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Fig. 8. Two-point velocity correlation function for the weaker 
absorbers, up to 1 000 kms -1 , in 100 kms -1 bins. The correla- 
tion function is shown by the solid line. For orientation £ = 
is marked with a dotted line. Dashed and dot-dashed lines rep- 
resent the lcr and 2<x Poisson errors, related to £ = 0, respec- 
tively. 
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Fig. 9. Two-point velocity correlation function for the weaker 
absorbers, up to 10 000 km s _1 , in 1000 kms -1 bins. Line sym- 
bols are the same as in Fig. [3] 

splittings, repeating this procedure 1000 times, and computing 
the average. 

When the clustering properties in multiple LOS are investi- 
gated, each quasar spectrum must first be analyzed separately. 
Subsequently, the individual numbers of pairs n b s and n s i m of 
the j analyzed LOS can be added in order to build a total £ 
based on improved statistics: 



^(Av)=^'"° b - (Ay) -l. 

£- =1 «sim,/(Av) 



(9) 



For varying intervals Az and AlogA^Hi, we counted n Q b s 
and simulated n^ m for every LOS. To investigate the cluster- 
ing on small scales, we chose Av max = 1 000 kms" 1 using a 
binning of 100 kms -1 , and we defined Av max = 10000 kms -1 
(in 1 000 kms -1 bins) for a study of the long-scale clustering 



Fig. 10. Two-point velocity correlation function for the 
stronger absorbers, up to 10000 kms -1 , in 1 000 kms -1 bins. 
Line symbols are the same as in Fig. [8] 
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Fig. 11. Two-point velocity correlation function for the 
stronger absorbers, up to 1 000 kms -1 , in 100 kms -1 bins. Line 
symbols are the same as in Fig. [8] 



properties. We calculated ^ according to © with j = 10 (be- 
cause the LOS to HS 0747+4259 was formally treated as two 
LOS: one for z < 1.44 and one for z > 1.56; compare with 
Table Qj. We estimated the error on the clustering signal with 
error propagation, assuming a Poisson error for n s \ m and setting 
cr„ obs = 0. We think that this is an acceptable approximation be- 
cause cr„ bs follows directly from crx or <x,, respectively, and the 
line positions are very well known compared with the values of 
the investigated line separations Az. Furthermore, « t, s and n s ; m 
have similar values in most cases. 

Since Cristiani et al. (1995) found that clustering is stronger 
for higher column density absorbers, we also looked for this 
trend. We took the conventional definition log A^hi > 13.64 
for the strong lines, and we chose log A^hi = 12.90 - 14.00 
for the low column density absorbers in accordance with the 
completeness limit of our column density distribution function 
(see Subsection l5.1> . 
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In Figs.[S]and[5]we show £(Av) for the low column density 
absorbers on scales up to 1 000 kms _I and 10000 kms -1 , re- 
spectively. We detect a marginal signal above the 2cr level only 
for Av = 150 kms -1 , 1500 kms~' and 6 500 kms -1 . Apart 
from that, £(Av) is always below the lcr level. The excess 
at Av = 150 kms" 1 = 0.32, lcr = 0.13) suggests a clus- 
tering of the weak lines over short distances. The signal at 
Av = 50 kms^ 1 = 0.02, lcr = 0.10) might be too weak be- 
cause some weaker lines remain unidentified due to blends: re- 
solving these blends into superpositions of several lines would 
result in far more line pairs at low Av. 

As expected, we also found slight clustering signals for the 
high column density absorbers, however only over short dis- 
tances: £ = 1.04 > 2cr for Av = 50 kms" 1 and £ = 0.36 > lcr 
for Av = 150 kms -1 (Fig. II li . At higher distances (up to 
1 000 kms ), the strong absorbers are obviously uncorrelated. 
The large-scale investigation of the stronger lines provides 
marginal signals exceeding the lcr level at Av = 500 kms -1 
and A7 = 7 50010ns- 1 (Fig.[l0j. 

The clustering of the high and low column density ab- 
sorb ers on short scales co ntradicts our own results presented 
in Ijanknecht et alJ d2002l) . where no clustering signal at all 
was found in the LOS to HE 0515-4414 alone, independent of 
the investigated column density ranges and velocity splittings. 
Individual LOS might not offer a sufficient statistical basis to 
detect a weak clustering signal. This is supported by Fig. [21 
where we plot the correlation function £ for some individual 
LOS of our sample. Although the four LOS with the highest 
numbers of strong lines are considered (> 40), occasionally 
an at most marginal excess above the lcr level can be seen for 
Av < 200 km s _1 , but no value at all above the 2cr level. 

Most studies found significant clustering on short scales, 
even to some extent for the low column density fines 
JKim et allboOlUCristiani et allll995llUlmerlll996llHu et all 
1199.4 . 

Previous observations sug gested a clustering signal increas - 
ing with decreasing redshift llKim et all 1200 ll llJlmer. 1996). 
This is made plausible by the fact that the structures in the uni- 
verse develop by gravitational pull over the course of time. In 
addition, at low reds hift > 1 /3 of all Ly a absorbers are asso- 
ciated with galaxies dLanzetta et allll995l) which are known to 
cluster on these scales. 

We also searched for a dependence of the correlation func- 
tion on z. For this purpose, we defined three z intervals of 
the same width (0.5 - 1.0, 1.0 - 1.5 and 1.5 - 2.0). As de- 
scribed above we computed £ and estimated the lcr and 2cr 
errors separately for each z interval and for distances up to 
Avmax = 1000 km s -1 in 100 kms~' bins. The results are il- 
lustrated in Figs. and [21 For the weaker lines, the £ « 
value for Av < 100 kms -1 is again probably caused by blend- 
ing effects, while for Av > 200 kms -1 no clustering signal 
can be recognized. For Av = 150 kms -1 the strongest cluster- 
ing is found in the highest redshift range: £ = 0.44 > 2cr. At 
lower z, £ is roughly at the lcr level (£ = 0.18, lcr = 0.19 for 
z = 1.0- 1.5; £ = 0.33, lcr = 0.30 for z = 0.5 - 1.0), i.e., for 
z < 1.5, £ doesn't change relative to cr. 

The high column density lines show a similar picture 
for Av = 50 kms -1 : starting at a ~ 2cr level (£ = 1.40, 
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Fig. 12. Two-point velocity correlation function for the 
stronger absorbers for individual quasar LOS, in 100 kms 1 
bins. Line symbols are the same as in Fig. [8] When there are 
no pairs in a specific bin, formally « Q b s = 0, £ = -1, and the 
1 cr and 2 cr Poisson errors are 0, following error propagation 
as was done for all other bins. 



lcr = 0.79), the clustering degree decreases strongly with de- 
creasing redshift. However, the opposite trend can be seen at 
Av = 150 kms -1 where the signal increases from ^ « to 
£ a 1 .5 cr when going from the highest to the middle z interval. 
Besides, neither the high nor the low column density absorbers 
show a clustering signal lying clearly above the 2cr level. 

In summary, no uniform evolutionary behaviour of £ can 
be established. The rising correlation of the absorbers with 
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Fig. 13. Evolution of the two-point velocity correlation func- Fig. 14. Evolution of the two-point velocity correlation func- 
tion for the weaker absorbers, in 100 kms bins, illustrated tion for the stronger absorbers, in 100 kms -1 bins, illustrated 
with three different z intervals. Line symbols are the same as in with three different z intervals. Line symbols are the same as in 



Fig. El 



Fig.E] 



decre asing redshift suggested bv lKim et al.l d200ll) and lUlmeJ 
( 1996) cannot be confirmed. 

5.4. Evolution of the number density 

As for the clustering properties, a different behaviour is also 
expected for the number density evolution of the stronger and 
weaker Ly a lines. We analyze this evolution using the two 
column density ranges from Subsection l5.3l separatelv. 

In practice, the value for the differential number density ^ 
can be deduced by choosing an interval width Az and count- 
ing the number of absorption lines n in the individual intervals. 
When the line samples of several LOS which in general cover 
different z regions are combined, it must be taken into account, 
of course, that a varying number of LOS «l falls into a cer- 
tain redshift interval. (Here, «l can be estimated from Fig.Q. 
Thus, we defined a density 4? per LOS. Since the fluctuations 



between the LOS are large, only z intervals with n L > 1 were 
considered. This is the case for 0.7 < z < 1.9. 

We computed ^ for the two sub-groups as a function of z. 
On the basis of the Levenberg Marquard algorithm, we fit the 
power law {1} to the data points of ^(z) and determined 
and y in this way. 

Figs. H5landll6lshow the data points for the low and high 
column density absorbers, respectively. In addition to the data 
points, the best fits to them with their 2cr confidence limits are 
presented. 

For the weak absorbers, the number density decreases with 
decreasing redshift with y = 0.7 '4 + 0.31. The outliers in the 
interval z = 0.7 - 0.8 might be due to the small size of our 
sample in this region (see Figs. H an d|3. The stronger lines 
evolve faster with an exponent y — 1 .50 + 0.45 which is only 
just consistent (within lcr) with the exponent of the weak lines. 

Fig. gives an extensive comparison of the results 
of this work for the weaker absorbers (here defined by 
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Fig. 15. Number density evolution of the weak Ly a absorbers 
(12.90 < log A H i < 14.00). Given are the data points (binned 
with Az = 0.05), the best fit and the 2<x confidence limits. 
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Fig. 16. Same as in Fig. ^] for the strong Ly a absorbers 
(log N m > 13.64). 



logA^Hi = 13.10- 14.00 to facilitate comparison with other 
results) with different contributions to the number density 
evolution from the literature measured in varying z regions. 
Fig-EDdoes the same for the stronger absorbers. In the double- 
logarithmic presentations the symbols correspond to the num- 
ber densities measured in different studies. The Ly a lines of 
this work are binned in intervals of Az = 0.2 in order to have in- 
terval sizes comparable with the other investigations. The data 
points of this study are shown with filled circles; the error bars 
result from a 1 <x poisson error of n. 

Regarding the evolution of the weaker absorbers (Fig.ll7L 
it becomes clear that the LOS analyzed here fall into a z region 
which has not yet been examined so far, due to a lack of high- 
resolution UV data. The best fit to the data points of this work 
gives j = 0.78 + 0.27. 

The decrease of the line density with decreasing redshift 
is obviously decelerated in our investigated z interval, as is re- 
vealed from a co mparison with higher redshifts: for z > 1.5, 
iKim et all J2002I) derived y = 1.18 + 0.14 from a best fit to 
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Fig. 17. Comparison of different studies of the number density 
evolution of the Ly a absorbers with 13.10 < log A^hi < 14.00. 
The Ly a lines of this work are binned in Az = 0.2 intervals. 
The solid line represents the best fit to the data points of this 
work, while the dotted line is an extrapolation of the fit for 
z — » 0. The dashed line shows the bes t fit to the data points in 
the higher redshift range, taken from lKim et alJ d2002l) . 
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Fig. 18. Comparison of different studies of the number den- 
sity evolution of the Ly a absorbers with log Ahi > 13.64. 
Details as in Fig. ^] The additionally given dashed line for 
log ( 1 + z) = - 0.4 repre sents the best fit to the data points 
from Wevmann et al .(1998). The data points from Penton et al. 
(2004), Savaglio et al. (1999) and Williger et al. (2003) apply 
to log A H i > 14.00. 



the data points of different studies (presented as dashed line in 

Fig-EJ- Excluding two LOS, they found y = 1.42+0.16 which 

is no longer in agreement with the result found here. 

Interestingly, an extrapolation of our fit line for z — > is 

consistent within lcr with the study of Pent on et al ] (l2004l) for 

the local universe (see dotted line in Fig. ll7l >. wh ile an extrap- 
I ' I j i ° — I I 

olation of the Kim et al. (2002) fit line misses the Penton et al. 

( 2004 ) data point by w 2cr. dWillieer et alJ J2003h examine only 
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a single LOS). We interpret this as further evidence for a chang- 
ing evolution of the lower column density lines in the interval 
Z = 1 - 2. 

Furthermore, the comparison with the investigations at 
higher and at lower z suggests a continuous and moderate tran- 
sition to the present, rather than a sharp break. 

As already derived from Figs. ^] and ^] the high 
column density absorbers show a steeper gradient in the evolu- 
tion diagram with y — 1.66 + 0.57 (Fig.ll8>. A remarkable data 
point occupies the interval z — 0.9 - 1 . 1 . Without this outlier the 
evolution index would be much higher and its error much lower 
(y = 2.11 + 0.21). A detailed analysis of this redshift range 
which is covered by six of the nine LOS (see Fig.^ shows that 
the high line density is not caused by a fluctuation in an indi- 
vidual LOS, but that the ^ lies clearly above the level of the 
adjacent intervals in nearly all six LOS. Obviously we detect a 
coincidental line accumulation within this interval. 

Just as for the weaker lines, our y is lower compare d with 
the results of investigations at higher z. Ki m et alJll2002l) found 
y = 2.47 + 0.18 by fitting a straight line to the data points of 
varying studies for z > 1.5 which is inconsistent within lcr 
with the value computed here. In addition, the number density 
continuously decreases in the analyzed range z = 0.7 - 1.9 
(apart from the outlier) without showing the freque ntly postu- 
lated break in the evolution at z = 1.5 - 1.7 (|Wevmann et all 



ll998UDobrzvcki et alll2002Hlmpev et ailll996l) . Possibly this 
slow-down in the evol ution, which was also successfully sim - 
ulated by theoreticians jTheuns et ailll998tlDave et all fl999). 
takes place in a later phase (z < 0.7) beyond the redshift region 
examined here. An indication for this might be that at lower z 
the thinning out of the absorbers obviously continues to decel- 
erate: in contrast to the weaker absorbers, an extrapolation of 
our fit line for z — > would lead to a distinctly too low ab- 
sorber density a t z = , compared with Pen ton et alJ J2004I) and 
IWevmann et aD J 19981) (see Fig.M^, 

If we fit the data points of Wevman n et alJ {l298) with a 
straight line (shown in Fig.^J y = 0.13 + 0.06), it does not in- 
tersect our fit line before z = 0.6. Also, leaving out our outlier 
leads to a point of intersection at z — 0.8. Thus, in the case that 
the break in the number density evolution took place, the red- 
shift value at which it occurred ha s to be corrected downw ards. 

However, it must be noted that Wevma nn et alJ dl998 ) used 
the equivalent width Wa instead of the column density as the 
observable. The conversion formula 

Merest [A] = 8.85 • 10- 21 (a ks1 [A]) 2 / N [cnT 2 ] (10) 

(/ oscillator strength) only applies to unsaturated absorption 
lines and depends on the Doppler parameter. Occasionally, 
b — 25 km s" 1 is assumed in the literature dPenton et all 
l2000l 12004 . Then the l ower limit for the equi valent width 
W^rat = 240 mA used bv lWevmann et alJ Jl998h corresponds 
to a minimum column density log A^Hi.min = 14.00. Applying 
this higher cut-off to the evolution diagram of the stronger 
absorbers would shift our data points in Fig. 1181 significantly 
downwards. As a re sult of this, the i ntersection of our fit line 
with the one to the IWevmann et alJ dl998l) data points (and 
therefore the evolution break estimated from it) had to be cor- 
rected to a higher redshift. For this reason, the uncertainty in the 



Wa-Nui conversion makes the exact localization of the slow- 
down difficult to determine. 

Ne vertheless, newer studies of lKim et alJfcOO lLlKim et alJ 
ll2002l) (whose fit line hits that of th e IWevmann et alJ jT 998) 
data points at z = 1.1; see Fig. II 81 or Penton et al. (2004) also 
favoured the break in the evolution taking place at no higher 
redshift than z - 1.0-1.2. Supported by these results, we con- 
clude that a transition in the evolution of the higher column 
density lines indeed took place, but much later than assumed 
before. In simulations the sudden slow-down in the evolution 
results directly from the decline of the UV background at z — 2 
jDaveetallll999h. However, the QSO dominated UV back- 
ground of lHaardt & Madaul i 19961) is usually assumed in these 
calculations, and thus the contribution of the galaxies is possi- 
bly underestimated. With the latter as main source of the UV 
radiation field for z < 2, the radiation field would decrease 
more slowly and the break would naturally be corrected to a 
lower z. 



All results for the evolution of the number density of Ly a 
absorbers in this subsection have to be considered in view of 
the high scatter of the data points. This scatter could not be sig- 
nificantly reduced compared with our stu dy of the single LOS 
to HE 05 15-4414 Jjanknecht et alll2002l) . though we enhanced 
the examined total redshift path length Az by roughly a factor 6. 
For all nine observed QSOs we obtained high-resolution spec- 
tra (R > 30 000) where the Ly a forest is completely resolved 
(because^Hi > 10 km s _I ). Therefore, we suppose that the vari- 
ations from LOS to LOS are caused physically (rather than me- 
thodically or statistically), reflecting a s trongly inhomoge nous 
appearance of the Ly a forest. Indeed, iKim et al.l ( I2002J) and 
iTvtler et al.l J2004ll (for z < 2.5 and for z = 1.9, respectively) 
and Imnev et al. (1999) for the local universe (z < 0.2) found 
hints for a strong cosmic variance too. 

The reason for the large variations of ^ between individual 
LOS falling into the same redshift ranges is probably the ad- 
vanced stage of structure formation in the universe at z = 2. The 
structures are gravitationally evolved: over- and underdensities 
have been formed. In individual LOS this can be recognized 
in the form of strongly varying number densities, in particular 
of the high column density absorbers. The structure formation 
obviously starts to dominate over the Hubble expansion, which 
causes the general decline of the number density. 



5.5. Effective optical depth atz~2 

We examined the effective optical depth r e ff and the average 
normalized flux (F norm ), which are related to each other by the 
equation 



Teff = -In 



= -In (F norm > , 



(ID 



where F cont is the continuum flux. We calculated r e fr 
and (Fnorm) for the maximum redshift range of our LOS, 
2=1.81-1.91. The LOS to three of our quasars fall into the 
chosen interval: those to HE 2225-2258, HE 0429-4901 and 
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HS 0747+4259 (see Fig.[T}. Therefore, we used the flux distri- 
butions of the corresponding spectral parts for the calculation 
of (F norm ). Because we were interested in the Hi absorption of 
the low density IGM, we defined wavelength regions around 
the metal lines occuring in the regarded z interval and elimi- 
nated these regions from the three spectra. The remaining flux 
distribution F(A) was divided by the respective continuum flux 
^cont(^) which is an output parameter of CANDALF. For the 
resulting normalized flux averaged over the wavelength posi- 
tions of all three spectral parts we derived 

<F nonn > = 0.895 ± 0.020 (12) 

for < z >~ 1.86, i.e. 10.5% of the radiative flux is absorbed 
by the diffuse intergalactic medium at this redshift. The error 
of (^norm) follows from Tytler et al. (2004), eq. (30): since we 
have only three LOS, the total error of our average normalized 
flux is completely dominated by the large scale structure and 
will probably be > 0.02. With (1111 and (II 21 we compute 

r eff = 0.111 +0.020 (13) 

for the effective op t ical d epth. 

iKirkman et al.l d2005l) investigated the Hi opacity of the 
IGM at 1.6 < z < 3.2. At z=1.86, the scaling in their Fig. 5 
gave (F aorm ) = 0.889 with an error of roughly 0.01. Thus, our 
result agrees with theirs. 

6. Conclusions 

We have analyzed the combined sample of all Ly a absorp- 
tion lines which were detected in the spectra of nine quasars 
over a total redshift path length Az - 5.176 within the 
range z = 0.5 - 1.9. We classified 1325 absorption features as 
Ly a lines. Examining the distributions of their fit parameters 
log Njji, b and z, we find the following: 

1 . The column density distribution might be complete down 
to log Nm = 12.90. It can be approximated over roughly 
three orders of magnitude (log Nui = 12.90 - 15.70) by a 
simple power law. The slope /? = 1.60 + 0.03 is consistent 
with other studies in comparable redshift intervals. We do 
not see any change in /3 with z. 

2. The distribution of the Doppler widths has the expected 
form of nearly Gaussian with an additional expanded tail 
to high b and an average value b = (34 + 22)kms~ 1 . We 
cannot detect any evolution of the Doppler parameter dis- 
tribution with z in the examined redshift phase. 

3. The weaker (log JV H i = 12.90 - 14.00) as well as the 
stronger Ly a lines (log A^hi > 13.64) show marginal 
clustering with a 2cr significance on short velocity inter- 
vals (Av < 200 kms~' and Av < 100 kms~', respectively). 
However, the clustering signal is too weak to be seen in 
individual LOS. We do not see any change with column 
density or redshift. 

4. Using the customary power law 4s oc (1 +z) y , the evolution 
of the number density of the weaker Ly a lines declines 
with y - 0.74 + 0.31 with decreasing z within the range 
z - 0.1 - 1.9. From comparisons with investigations of 



other authors at higher and lower redshifts we conclude that 
the decrease of the weaker lines is decelerated in the phase 
Z = 1 — 2, turning into a nearly flat evolution for z — » 0, 
without showing any hint for a sharp break in the evolu- 
tion. 

5. The number density of the stronger absorbers also de- 
creases with decreasing z (y = 1.50+0.45 for z = 0.7-1.9), 
faster than for the weaker lines, though the trends are con- 
sistent with each other within Icr. 

6. The decline of the line density of the stronger ab- 
sorbers is also decelerated compared with higher redshifts 
(z — 1.5 - 4.0). The break in the evolution predicted for 
z = 1.5-1.7 cannot be detected, however. The deduced 
z dependence of the number density as well as comparisons 
with other studies for the local universe rather suggest a 
later slow-down (z < 0.7), followed by an approximately 
flat evolution. The precise z value of the break in the evolu- 
tion probably has to be corrected somewhat upwards due to 
the uncertain conversio n between equivalent width (used by 
IWevmann et alJ Jl998h as observable) and column density. 

7. Though having analyzed nine high-resolution (R > 30 000) 
quasar spectra altogether, the scatter of the data points in 
the diagrams of the number density evolution is large. As 
Tytler et al. (2004) showed, the strong variations between 
the LOS are consistent with the effects of normal structure 
formation. 

8. An analysis of the flux distribution of the three quasars with 
LOS in the upper redshift range (z = 1.81 - 1.91) gives 
(^norm) = 0.895 + 0.020 for the average normalized flux 
and T e ff = 0.1 1 1 + 0.020 for the effective optical depth, in 
good agreement with the values from Kirkman et al. (2004). 
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